1 


Jn?:  m  i 


AD-A206  683  3ocumentatiow 


PAGE 


form  Approved 
OMiNo.  0704^1  U 


tb  RESTRiaiVE  MARKINCS 


2a.  security  classification  authori 


2b.  DECLASSIFICATION  /  OOWNC 


I  'll 


4.  PERFORMING  ORGANIZATION  R 


6a.  NAME  OF  PERFORMING  ORGANIZATION 

Princeton  University 


6c  ADDRESS  (City,  Stata,  and  Zf^Codtj 

Princeton  New  Jersey  08544 


3  DISTRIBUTION /AVAILA6IUTY  OF  REPORT 

Approved  for  public  release; 
distribution  is  unlimited. 


S.  MONITORING  organization  REPORT  NUMIER(S) 

AhObK-lli-  D  -  0  4  j  r 


7a.  NAME  OF  MONITORING  ORGANIZATION 

AFOSR/NA 


7b.  ADDRESS  (Oty.  Statt,  and  ZIP  Codb) 

Building  410/  Bolling  AFB  DC 
20332-6448 


8a.  NAME  OF  FUNDING /SPONSORING 
ORGANIZATION 

AFOSR/NA 


•c  ADDRESS  (Ofy.  Sta«a.  and  ZIP  Codb) 
Building  410/  Bolling  AFB  DC 
20332-6448 


1 1 .  TITLE  (indudp  Soturtty  Ostdfkation) 

Lumped  Model  Generation  and  EvaluationrSensitivity  and  Lie  Algebraic  Techniques 
'  With  Applications  to  Combustion  _ 


12.  PERSONAL  AUTM0R(S) 

H.  Rabitz.  F.  L.  Dryer  and  R.  Yetter  _ _  . 


tia.  TYPE  OF  REPORT  1 3b.  TIME  COVERED  114.  DATE  OF  REPORT  (Ytpr, Month, Opy)  15.  PAGE  COUNT 

Final  Technical  FROM  9/10/85  TO  11/30/8^  89,03,03 _ 9 


16.  SUPPLEMENTARY  NOTATION 


COSATI  COOES 


GROUP  SUB-GROUP 


IB.  SUBJECT  TERMS  (CondniM  on  rtvora*  if  notoaary  and  idanmy  by  btocR  numbar) 
coobustion  modelling,  senslclvlcy  analysis, 
system  lumping  and  reduction.  Lie  groups 


19.  abstract  (Cendmis  on  rovom  if  natoaary  and 


by  block  rwmbar) 


This  program  dealt  with  the  development  and  application  of  new  approaches  for  producing 
and  evaluating  seml-emplrlcal  (lumped  parameter)  models  of  physical  processes.  Pro¬ 
cedures  using  local  sensitivity  gradient  methods  were  used  to  study  the  existing  lumped 
kinetic  models  for  the  moist  carbon  monoxide  oxidation  to  show  that  transport  processes 
can  cause  oversimplified  lumped  models  derived  from  homogeneous  kinetics  to  fail  when 
applied  to  flame  propagation  systems.  New  models  are  under  developioent  which  will 
Include  the  appropriate  level  of  detail.  A  Lie  group  formalism  was  developed  to 
address  global  parameter  space  mapping  Issues  for  first  order  differential  equations. 
The  rigorous  criteria  for  the  existence  of  exact  lumping  by  linear  projective  trans¬ 
formations  was  also  established,  ^ . 


20^  DISTRIBUTION /AVAIIABIUTY  OF  ABSTRACT 
S  UNCLASSIFIEO/UNLIMITEO  S  SAME  A$  RFT 


224.  NAME  OF  RESFONSIBLE  INDIVIDUAL 

Julian  M  Tishkoff 


OOForm1473/  JUN  88  ,  ^ 


TIC  USERS 


|22b  TELEPMONLf//K/«d4iWf4Codt)  1 2Ze^9FFiCE  y  VSC 

1  (202)767-0465  _ I  AFOSR/NA 


fravioui  aOitiont  art  obtoiato. 

^  : 


riSn!3EnrCT52!7S5H 


Unclasaif led 


1 

APOSR  TR.  8  9-0418 

Abstract 

Research  toider  this  program  is  aimed  at  setting  out  criteria  for  the  ability 
to  lump  or  reduce  complex  combustion  kinetics  systems.  Research  during  this  first 
phase  of  the  project  primarily  focused  on  well  stirred  chemically  reacting  systems 
described  by  ordinary  differential  equations.  The  complexity  of  this  problem 
necessitates  a  multi-faceted  approach  and  three  principal  avenues  of  pursuit  were 
considered.  First,  the  tools  of  sensitivity  analysis  were  employed  to  Indicate 
particular  portions  of  a  chemical  mechanism  of  importance  as  well  as  those  which 
may  be  eliminated.  In  carrying  forth  this  approach,  an  interesting  and  surprising 
level  of  scaling  and  self  similarity  was  observed  amongst  the  sensitivity 
coefficients  strongly  implying  an  ability  to  lump  and  reduce  the  models.  Secondly, 
the  rigorous  criteria  for  lumping  under  a  class  of  linear  projective  trans¬ 
formations  were  established.  These  criteria  have  utilities  for  a  wide  variety  of 
nonlinear  kinetic  systems  as  well  as  other  applications.  Finally,  a  fully 
nonlinear  approach  to  lumping  and  general  system  reduction  was  developed  based  on 
Lie  group  techniques.  Each  of  these  three  approaches  has  its  own  realms  of 
application  and  the  research  under  this  effort  has  clearly  set  the  foundation  for 


further  development  aimed  ultimately  at  practical  applications. 
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Siinmiarv  of  Research  Activities 

Given  below  is  a  synopsis  of  the  research  carried  out  iinder  this  program. 

Each  section  summarizes  a  component  of  the  research.  The  superscripts  appearing  in 
the  section  titles  refer  to  specific  papers  which  have  appeared  on  these  subjects. 
These  references  are  summarized  on  page  8  of  this  report. 

1.  Scaling  Relations  and  Self-Similarity  Conditions  in  Strongly  Coupled 

Dynamical  Systems^ 

Dynamical  equations  arising  in  a  number  of  physical  areas  typically  involve 
many  dependent  variables  as  well  as  numerous  parameters.  In  some  cases  these 
equations  may  contain  a  dominant  dependent  variable  (e.g.,  the  temperature  in  flame 
systems,  etc.),  and  the  consequences  of  such  an  identification  are  examined  in  this 
paper.  Particular  emphasis  is  placed  on  the  behavior  of  dynamical  Green's 
functions  and  system  sensitivity  coefficients  with  respect  to  the  parameters 
residing  in  the  particular  model.  In  the  case  of  a  single  dominant  dependent 
variable  for  a  system  of  ordinary  differential  equations,  it  is  possible  to  reduce 
the  Green's  function  matrix  to  knowledge  of  one  of  its  columns  and  often  to  one 
independent  element.  Furthermore,  idien  there  are  N  dependent  variables  and  M 
parameters,  the  N  x  N  matrix  of  sensitivity  coefficients  reduces  to  knowledge  of 
only  two  characteristic  vectors  of  lengths  N  and  N,  respectively.  These  various 
reductions  are  referred  to  as  scaling  relations  and  self* similarity  conditions. 

The  consequences  of  a  dominant  dependent  variable  are  illustrated  with  examples 
drawn  from  various  areas  of  combustion  and  kinetics  modelling.  A  brief  discussion 
in  the  Appendix  is  also  presented  on  similar  organizing  principles  in  multidominant 
dependent  variable  systems  and  cases  described  by  partial  differential  equations. 
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2.  A  General  Analysis  of  Exact  Lumping  in  Chemical  Kinetics^ 

A  general  analysis  of  exact  Ituaping  is  presented.  This  analysis  can  be 
applied  to  any  reaction  system  with  n  species  described  by  a  set  of  first  order 
differential  equations  dy/dt  -  f(y).  where  y  is  an  n-dimensional  vector;  f(y)  is  an 
arbitrary  n-dimensional  function  vector.  Here  we  consider  lumping  by  means  of  an 

nxn  real  constant  matrix  M  with  rank  n(n<n) .  It  is  found  that  a  reaction  system  is 

* 

exactly  lumpable  if  and  only  if  there  exist  nontrivial  fixed  invariant  subspaces  H 
of  the  transpose  of  the  Jacobian  matrix  Jt (y)  of  f(y),  no  matter  what  value  y 
takes  and  the  corresponding  eigenvalues  are  the  same  for  Jt (y)  and  Jt (M  My) .  Here 
the  rows  of  M  are  the  basis  vectors  of  M  and  M  is  any  generalized  inverse  of  M 
satisfying  M  M  -  !„  with  being  the  n-identity  matrix.  The  fixed  invariant 
subspaces  of  J^ (y)  can  be  obtained  either  from  the  simultaneously  invariant 
subspaces  of  all  A]^,  where  the  A^^'s  form  the  basis  of  decomposition  of  J^  (y) ,  or  by 
determining  the  fixed  KerCw^CJ^  (y)  -  ^  ^n  "  ^aj  (y) 

+  (J^(y))^]  )  I  where  aj  ±  irj  are  the  real  and  nonreal  eigenvalues  of  J^(y) 

and  A£,(7j  and  rj  are  usually  functions  of  yj^  and  rj  are  non-negative  integers.  The 
kinetic  equations  of  the  lumped  system  can  be  described  as  dy/dt  -  Mf(My).  This 
method  is  illustrated  by  some  simple  examples. 

3.  Identifiability  and  Distinguishability  of  First-order  Reaction  Systems^ 

By  following  the  kinetics  of  a  reaction  through  the  use  of  certain  classes  of 
measurable  quantities  instead  of  the  concentration  of  all  species,  neither  the 
parameter  values  nor  the  reaction  scheme  is  necessarily  iinique.  Identifiability 
deals  with  the  problem  of  determining  idiether  an  experiment  is  able  to  supply  the 
desired  information  on  the  parameters  of  an  assumed  kinetic  model,  whereas 
indistinguishability  means  that  two  different  reaction  schemes  generate  the  same 
values  for  the  observed  quantities  in  any  possible  experiment.  This  paper  examines 
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these  Issues  for  the  case  of  first-order  reaction  systems,  and  both  problems  are 
solved  by  the  same  analytical  tools.  The  method  Involving  Laplace  transforms  Is 
conceptually  simple,  is  easy  to  apply,  and  is  also  used  to  derive  simple  rules  to 
test  dlstingulshablllty  of  reaction  schemes.  Another  approach  based  on  similarity 
transformations  Is  used  to  generate  all  the  first -order  reaction  schemes  that  are 
indistinguishable  from  a  given  one. 

4.  Parameter  Space  Mapping  of  First-Order  Linear  Ordinary 

Differential  Equations^ 

A  Lie  group  formalism  for  global  parameter  space  mapping  of  ordinary 
differential  equations,  described  In  a  preceding  paper.  Is  developed  here.  The 
need  for  such  mapping  arises  in  a  variety  of  physical  contexts.  The  procedure  is 
demonstrated  on,  but  not  restricted  to,  the  system  of  coupled  equations  x  -  cx. 

The  Lie  group  generators  are  obtained  exactly  from  both  the  time -Independent  and 
time -dependent  generating  equations.  The  transformations  obtained  from  these 
generators  leave  the  system  of  differential  equations  invariant.  The 
tlme-lndpendent  transformations  map  any  solution  of  the  linear  system  x  -  xc  Into 
any  other  solution  with  the  same  frequency  or  time  constant.  The  time -dependent 
transformations  interconvert  solutions  with  different  frequencies  and/or  small  time 
constants.  Any  solution  of  x  -  xc  can  be  mapped  Into  any  solution  of  x  -  c'x. 

Thus  the  behavior  of  x(t)  can  be  examined  as  a  function  of  changes  In  any  of  the 
system  parameters,  or  in  the  Initial  conditions,  x(0).  As  an  example,  one  of  the 
time -dependent  mappings  demonstrates  the  continuous  transformation  of  oscillator 
solutions  into  nonosclllatory  ones,  and  vice  versa. 
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i.  Global  Sensitivity  Analysis  of  Nonlinear  Chemical  Kinetics  Equations 

Using  Lie  Groups:  I.  Deterntnatlon  of  One -parameter  Groups^ 

We  introduce  one -parameter  groups  of  transformations  that  affect  wide-ranging 
changes  in  the  rate  constants  and  Input/output  fluxes  of  homogeneoiis  chemical 
reactions  involving  an  arbitrary  number  of  species  in  reactions  of  zero,  first  and 
second  order. 

Each  one -parameter  group  is  required  to  convert  every  solution  of  such 
elementary  rate  equations  into  corresponding  solutions  of  a  one -parameter  family  of 
altered  elementary  rate  equations. 

The  generators  of  all  allowed  one-parameter  groups  are  obtained  for  systems 
with  N  species  using  an  algorithm  which  exactly  determines  their  action  on  the  rate 
constants,  and  either  exactly  determines  or  systematically  approximates  their 
action  on  the  concentration.  Compounding  the  one-parameter  groups  yields  all 
many-parameter  groups  of  smooth  time -independent  transformations  that  interconvert 
elementary  rate  equations  and  their  solutions. 

6.  Complications  of  One-Step  Kinetics  for  Moist  CO  Oxidation* 

The  one- step  reaction  mechanism,  CO  +  Oj  -►  COj  with 

dICO]/dt  -  -kov[C0]®[H20]'’[02]® 

which  is  frequently  used  in  combustion  problems  when  simplified  chemistry  is 
necessary,  is  numerically  studied  in  order  to  (i)  define  its  limitations  (and 
therefore  usage)  and  (ii)  understand  the  chemical  and  physical  reasons  for  these 
limitations.  The  analysis  is  carried  out  with  the  aid  of  a  validated 
contprehenslve,  elementary  reaction  mechanism  for  moist  CO  oxidation  and  by 
specialized  sensitivity  coefficients  which  correlate  the  parameters  of  the  global 
model  to  the  parameters  of  the  elementary  model.  The  results  confirm  many  of  the 
previous,  empirically  derived,  literature  models  and  show  the  overall  rate 
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constant,  as  a  function  of  temperature,  to  exhibit  non-Arrhenius  kinetics  and  to  be 
dependent  on  pressure  and  mixture  equivalence  ratio.  More  Importantly,  models 
derived  from  temporally  reacting  systems  are  shown  to  be  Improper  for  txse  in 
modelling  systems  reacting  with  transport  phenomena.  The  specialized  sensitivity 
coefficients  are  tised  to  explain  these  complex  behaviors  In  the  overall  model.  For 
the  temporal  system,  these  coefficients  show  that  the  global  model  must  be  able  to 
account  for  dissociation  and  equilibration  at  high  temperatures,  for  explosion 
phenomena  in  the  Intermediate  temperatures,  and  for  reaction  of  carbon  monoxide 
with  both  the  hydroxyl  radical  and  hydroperoxy  radical  at  low  temperatures. 

Lastly,  methods  for  modifying  the  existing  model  or  for  developing  a  new  model  are 
suggested. 


7 .  Elementary  Model  Reduction  of  Carbon  Monoxide  Kinetics 

In  our  previous  studies,  normalized  elementary  sensitivity  coefficients  have 
been  used  to  determine  the  relative  Importance  of  elementary  reactions  or  certain 
groups  of  reactions  In  comprehensive  mechanisms.  In  this  study,  we  have  extended 
this  methodology  by  using  the  principal  component  analysis  method  In  order  to 
systematically  reduce  the  size  of  the  original  comprehensive  mechanism.  Briefly, 
this  methodology  Is  based  on  a  least  squares  fit  approach  by  first  defining  the 
response  function,  Q,  as: 


Q(a) 


t  “  r  Ci(tj,a)  -  CiCtj.o*) 
J-1  iil  I-  Ci(tj.2*) 


where  a*  Is  the  nominal  values  of  the  parameters,  then  by  introducing  the  classical 
Gauss -approximation  to  yield  Q  (q)  Q(o)  -  (Aa)^s'^S(Aa)  where  aj  - 
finally  by  performing  an  eigenvalue -eigenvector  decomposition  on  the  resulting 
cross-product  matrix  S^S.  Eigenvectors  corresponding  to  small  eigenvalues 
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Indicate  unimportant  reactions,  thereby  enabling  one  to  optimally  reduce  the 
mechanism. 

Along  these  lines,  we  have  continued  our  previous  work  by  applying  this 
methodology  to  the  CO/H2/O2  reaction  mechanism.  A  large  number  of  Isothermal 
temporal  problems  were  ntunerlcally  run  to  generate  a  data  base  which  would  be 
representative  of  combustion  environment:^.  The  data  base  covered  a  temperature 
range  from  800  to  1800  K,  several  equivalence  ratios  from  lean  to  rich  conditions, 
and  several  pressures.  The  principal  component  analysis  was  applied  to  this  data 
base  to  determine  the  minimum  reaction  set  that  would  reproduce  all  the  original 
species  concentrations  within  2%. 

The  results  showed  that  the  original  52  reaction  mechanisms  could  be 
successfully  reduced  to  one  consisting  of  27  reactions  while  retaining  all  12 
species  In  the  model. 

Obviously,  this  reduction  is  still  not  practical  for  use  in  large 
multidimensional  codes.  The  necessary  further  reductions  are  proceeding  along 
several  directions.  First,  the  constraint  of  retaining  all  species  will  be 
lifted.  Our  earlier  research  has  shown  that  in  addition  to  the  major  reactants  and 
products,  two  Intermediate  species  are  necessary  In  the  model.  Secondly,  we  have 
also  found  that  In  more  complex  environments,  such  as  adiabatic  premixed  flames, 
the  underlying  chemical  processes  are  much  more  coupled  (namely  through  the  heat 
release  of  the  reaction)  and  hence,  such  problems  are  anticipated  to  be  more 
directly  Ivimpable. 
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